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4  INTRODUCTION 


In  the  United  States,  close  to  5  -  20  %  of  newly  diagnosed  breast  cancer  and  10  -  30%  of 
all  primary  breast  cancer  is  di  agnosed  as  locally  advanced  breast  cancer  (LABC)  [1,2]. 
Neoadjuvant  chem  otherapy  (system  ic  preope  rative  chemotherapy)  is  currently  the 
standard  of  care  for  LABC  [3,4].  When  compared  with  adjuvant  chem  otherapy 

(postoperative  therapy),  neoadjuvant  chemotherapy  yields  similar  results  for  both  overall 
survival  (70%  for  both)  and  disease-free  surv  ival  (53%  adjuvant,  55%  neoadjuvant)  [5]. 
Thus,  the  postponement  of  surgery  does  not  affect  the  outcome  of  the  treatment  [5,  6].  In 
addition,  neoadjuvant  chem  otherapy  offers  considerable  benefits  to  the  patient  as  the 
treatment  can  shrink  the  turn  or  and  even  in  som  e  cases  offer  complete  path  ologic 
response  [3,  7].  This  reduction  in  turn  or  si  ze  increases  the  possi  bility  of  breast 
conservation  [3,  5-7].  Maxim  izing  the  conser  vation  of  breast  tissue  can  be  of  great 
personal  im  portance  for  the  se  lf-esteem  and  quality  of  liv  ing  of  the  patient  [6]. 
Neoadjuvant  chemotherapy  can  also  offer  an  ear  ly  indication  of  the  patient’s  response  to 
chemotherapy.  Consequently,  m  onitoring  turn  or  response  to  neoadjuvant  therapy  gives 
the  possibility  of  adjusting  th  e  treatment  if  the  patient  is  responding  poorly  or  not  at  all 
resulting  in  substantial  advantages  for  the  patient  [3,  6].  This  project  aims  at  establishing 
noninvasive  m  onitoring  of  neoadjuvant  ch  emotherapy  in  the  breast  using  subhannonic 
aided  pressure  estimation  (SHAPE;  U.S.  Patent  6,302,845). 

Generally  interstitial  fluid  pressure  (IFP)  is  10-30  mmHg  higher  in  cancerous  tissue  than 
in  normal  tissue  although  values  of  up  to  60  m  mHg  have  been  recorded  [8,  9].  Sim  ilarly, 
IFP  in  breast  cancer  tumors  has  been  shown  to  be  higher  than  that  of  surrounding  breast 
tissue  [9].  This  increase  is  believed  to  be  due  to  vascularity,  fibrosis  and  difference  in  the 
interstitial  m  atrix  in  tumors  and  it  can  result  in  poor  transport  of  therapeutic  drugs  to 
tumors  [8].  Taghian  et  al.  used  a  wick-in-n  eedle  technique  to  m  onitor  the  IFP  of  br  east 
cancer  before  and  after  neoadjuv  ant  chem  otherapy  with  two  drugs  used  consecutively 
[10].  When  used  as  a  first  drug  Paclitaxel  decreased  the  IF  P  by  36%  (  p=0.02)  whereas 
with  Doxorubicin  as  a  first  drug  there  was  only  8%  reduction  (p=0.41).  As  this  was  a 
hypothesis-generating  study  they  did  not  show  any  outcome  related  to  the  relationship 
between  IFP  and  therapy  response  [10].  Howe  ver,  the  level  of  IFP  has  been  shown  to 
predict  disease  free  survival  for  cervix  cancer  (34%  disease  free  survival  (DFS)  if  IFP  > 
19  mmHg,  68%  DFS  if  IFP  <  19  mmHg  (p  =  0.002  ))  [11].  Thus,  the  le  vel  of  interstitial 
fluid  pressure  (IFP)  in  breast  cancer  turn  ors  could  potentially  be  used  to  m  onitor  the 
response  to  neoadjuvant  chemotherapy. 

Contrast  agents  have  b  een  used  for  two  d  ecades  to  im  prove  visualization  in  ultrasound 
(US)  imaging  as  they  enhance  the  difference  in  reflectivity  between  tissues  [12].  Because 
of  the  diffe  rence  in  compressibility  be  tween  the  m  edium  and  the  microbubble  any 
changes  in  pressure  induce  changes  in  the  size  of  the  m  icrobubble  [13].  This  in  turn 
affects  the  reflectivity  and  resonant  freque  ncy  of  the  bubble  [13,  14].  In  subharmonic 

imaging  (SHI)  pulses  are  transmitted  at  a  frequency  fo  and  the  echoes  are  received  at  half 
that  frequency  fo/2.  SHI  has  been  showed  to  be  a  f  easible  option  for  contrast  enhanced 
imaging  due  to  subharmonic  generation  by  c  ontrast  agents  and  lim  ited  subhannonic 
generation  in  tissues  [15].  Our  group  cam  e  up  with  a  novel  technique,  SHAPE,  utilizing 
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microbubbles  and  the  subhann  onic  amplitude  of  the  scattered  signal  [13].  W  e  showed 
that  there  is  a  linea  r  relationship  between  the  hydrostatic  pressure  and  the  subh  annonic 
amplitude.  We  propose  the  use  of  SHAP  E  to  m  onitor  treatm  ent  response  by 

noninvasively  m  easuring  the  IFP  in  breast  tumors.  This  offers  several  benefits  to  the 
patient.  As  opposed  to  the  wick-in-needle  technique  SHAPE  is  noninvasive  and  does  not 
inflict  pain.  Furthenn  ore,  it  allow  s  for  an  e  arly  indication  of  responders  vs.  non¬ 
responders  and  thereby  makes  adjustments  to  therapy  easier.  Moreover,  SHAPE  has  been 
shown  to  have  a  favorable  signal-to-noise  ratio  so  the  subharm  onic  am  plitude  is  not 
affected  by  background  noise  [13]. 

The  optimal  contrast  agent  and  acoustic  parameters  for  SHAPE  will  be  established  using 
in  vitro  pulse-echo  m  easurements.  The  SHAPE  algorithm  will  then  be  designed  and 
implemented  on  a  commercial,  state-of-the-art  US  scanner  for  in  vivo  IFP  measurements. 
A  sim  ilar  algorithm  has  already  been  set  up  for  cardiac  SHAPE  and  thus  only  a  few 
adjustments  need  to  be  m  ade  to  i  mplement  SHAPE  f  or  breast  turn  ors  making  this  very 
cost-effective.  The  in  vivo  experiments  will  be  twofold.  Fi  rst,  athymic,  nude,  female  rats 
will  be  implanted  with  SKBR3,  MCF-7  or  BT474  human  breast  cancer  cells  and  SHAPE 
used  to  measure  IFP  and  calibrated  by  comparing  the  SHAPE  results  to  IFP 
measurements  obtained  with  an  invasive,  in  tra-compartmental  p  ressure  m  onitor  as  the 
gold  standard.  After  calibration,  human  xenograft  breast  tumors  in  athymic,  nude,  female 
rats  will  be  used  to  evalua  te  the  ability  of  SHAPE  to  track  changes  in  IFP  by  studying 
before  and  after  administration  of  a  chemotherapy  agent  (paclitaxel). 

Our  group  has  proposed  that  SHAPE  and  contrast  enhanced  US  i  maging  can  be  used  to 
measure  the  IFP  in  LABC  turn  ors,  thus,  making  it  possible  to  noninvasively  m  onitor  the 
tumor  response  to  neoadjuvant  chemotherapy.  This  m  ethod  would  be  a  considerable 
improvement  from  the  wick-in-needle  technique  currently  used  for  IFP  measurements  in 
LABC  and  allow  for  individualized  treatments  options. 


5  BODY 

The  hypothesis  of  this  project  is  that  IFP  in  breast  tumors  can  be  measured  noninvasively 
using  SHAPE  and  contrast  enhanced  US  thus  improving  the  m  onitoring  of  neoadjuvant 
chemotherapy.  To  investigate  this  prospect,  in  vitro  pulse-echo  experim  ents  will  be 
conducted  to  investigate  this  prospect  and  find  the  optim  al  contrast  agent  for  SHAPE. 
These  re  suits  will  then  be  used  to  im  plement  SHAPE  on  a  comm  ercial  scanner.  The 
scanner  will  be  used  for  in  vivo  studies  on  201  rats  with  tu  mor  xenografts  in  order  to 
calibrate  and  evaluate  SHAPE’  s  ability  to  m  onitor  response  to  neoadjuvant 
chemotherapy.  The  specific  tasks  of  the  project  (as  presented  in  the  original  Statement  of 
Work)  can  be  found  in  Appendix  I. 

First  an  outline  of  the  m  ethods  applied  will  be  given  followed  by  a  presentation  of  the 
results  to  d  ate.  Fin  ally,  the  cone  lusions  a  nd  future  directions  of  the  research  will  be 
discussed. 
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5. 1  Methods 

In  vitro  experiments 

Last  year  a  pulse-echo  system  was  c  onstructed  to  test  different  types  of  contrast  agents 
for  use  with  high  frequency  SHAPE.  The  subharmonic  a  mplitude  at  different  static 
pressures  was  m  easured  using  a  sealed  wa  ter  tank  capab  le  of  w  ithstanding  pressure 
changes  over  200  mmHg.  Single  elem  ent  transducers  with  center  frequencies  of  4  to  12 
MHz  were  used  as  trans  mitter  and  receiver.  The  pressure  inside  the  ta  nk  was  monitored 
by  a  pressure  gauge  (OMEGA  Engineering,  Stamford,  CT).  However,  the  acoustic 
pressures  required  to  obtain  reasonable  subhannonic  signals  was  very  high  (0.9  —1.2 
MPa).  When  further  experiments  were  conducted  at  lower  acoustic  pressures  (0.5  to  0.8 
MPa)  and/or  a  higher  transm  it  frequencie  s  (10  -  12  M  Hz),  we  discovered  that  the 
attenuation  of  the  acoustic  wi  ndow  of  the  water  tank  m  arkedly  influenced  results  (i.e., 
reduced  or  eliminated  the  subharmonic  signal  components).  This  was  b  ecause  the  water 
tank  em  ployed  was  originally  designed  to  withstand  pressures  up  to  200  mmHg  (i.e., 
cardiac  pressures)  rather  than  the  smaller  pressures  encountered  in  breast  tumor  IFPs  (up 
to  50  mmHg). 

Hence  it  was  decided  to  repeat  the  static  p  ressure  measurements  using  a  m  ore  suitable 
water  tank.  Sm  all  (10  m  1)  OptiCell  chain  bers  (Nunc,  Rochester,  NY)  have  been  used 
successfully  to  investigate  in  teractions  between  contrast  bubbles,  ultrasound  and  cancer 
cells  [16]  and  these  were,  therefore,  selected  as  the  new  container  for  static  pressure  tests. 
The  OptiCell  was  submerged  into  a  larger  water  tank  and  0.2  ml/1  of  Dcfinity  (Lantheus 
Medical  Imaging,  N  Billerica,  MA)  injected.  RF  signals  were  acquired  with  a  Sonix  RP 
ultrasound  scanner  (Ultrasonix,  Richmond,  BC,  Canada)  using  a  high  frequency,  linear 
array  and  compared  to  an  invasive  (needle-in-wick)  pressure  monitor  (Stryker,  Berkshire, 
UK)  as  the  reference  standard. 

A  2.5  c  m  tissue  m  imicking  phanto  m  was  placed  between  the  probe  and  the  OptiCell 
chamber  to  simulate  tissue  attenuation.  Two  different  transmit  frequencies  of  6.7  and  10 
MHz  were  considered.  Contras  t  echoes  were  received  at  h  alf  the  trans  mit  frequencies 
(i.e.,  3.35  and  5.0  MHz).  The  acoustic  output  power  was  varied  from  0  to  -20  dB  for  a  0 
mmHg  hydrostatic  pressure  in  order  to  establish  the  optimal  sensitivity  for  SHAPE.  Then 
the  eh  amber  pressu  re  was  va  ried  f  rom  0  to  50  mmHg  to  sim  ulate  IF  P  in  turn  ors  and 
acoustic  pressure  v  aried  from  -4  to  -14dB.  A  fter  data  retrieval  the  amplitude  of  the 
subharmonic  signal  component  was  extracted  using  MATLAB  7.0.4  (Mathworks,  Natick, 
MA).  Three  m  easurements  were  a  cquired  at  each  setting  and  linea  r  regression  analysis 
used  to  determ  ine  the  relationship  betw  een  hydrostatic  pressu  re  and  change  in 
subharmonic  amplitude.  All  statistical  analys  es  were  cond  ucted  using  Stata  9.0  (Stata 
Corporation,  College  Station,  TX). 

Moreover,  a  novel,  sim  ulation  m  odel  of  the  dynam  ics  of  a  n  encapsulated  m  icrobubble 
contrast  agent,  developed  as  part  of  a  previous  DOD  s  upported  project  [17],  was 
modified  in  order  to  account  for  a  mbient  pressure  variations  and  different  shell 
parameters  to  establish  the  optim  al  c  ontrast  m  icrobubble  for  SHAPE.  A  nonlinear 
extension  of  the  origin  al  visco  elastic  m  odel  was  pursued  by  considering  a  quadratic 
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elasticity  model  where  the  interfacial  elasticity  vary  linearly  with  area  fraction  as  well  as 
an  exponent. 


In  vivo  experiments 

Our  group  has  worked  in  partnership  with  U1  trasonix  Medical  Corporation  to  im  plement 
SHAPE  for  cardiac  use  on  a  state-o  f-the-art  commercial  scanner  Sonix  RP  (Ultrasonix 
Medical  Corporation,  R  ichmond,  BC,  Canada)  with  a  phased  array  (P  A4-2).  Several 
experiments  have  been  carried  out  in  canin  es  to  investigate  cardiac  SHAPE  supported  by 
funding  from  the  AHA.  RF  data  from  thes  e  experim  ents  was  analyzed  off-line  using 
Matlab.  The  software  developed  and  optimized  by  our  group  for  this  analysis  is  not  site- 
specific  and  will  also  b  e  used  to  analyze  the  data  we  will  acquir  e  from  in  vivo  breast 
SHAPE  in  rats  as  part  of  this  project. 

Finally,  an  opportunity  to  provide  a  proof-of- concept  of  the  use  of  SHAPE  for  estimating 
IFP  and  test  the  invasive  (needle  based)  Stryker  pressure  m  onitoring  system  (the 
reference  standard)  presented  itself.  As  part  of  an  ongoing  NIH  study,  a  unique,  naturally 
occurring  tu  mor  m  odel,  the  Sine  lair  swin  e  with  melano  ma,  was  being  studied.  We 
obtained  IF  P  m  easurements  from  the  turn  ors  and  surrounding  tissue  using  the  S  tryker 
needle  based  system ,  which  provided  us  with  a  chance  to  as  sess  the  d  ependence  of  this 
technique  on  the  angle  between  the  needle  and  the  tissue.  Moreover,  subharm  onic 
signals  were  acquired  d  uring  an  infusion  of  Definity  (7.5  ml/l/min)  with  the  Sonix  RP 
and  a  linear  array.  Data  were  obtained  at  6.7  and  10  MHz  (i.e.,  subharmonic  frequencies 
of  3.75  and  5.0  MHz,  respectively)  with  acoustic  outputs  of  -4  and  -8  dB.  Five  (5)  swine 
were  studied  (one  melanoma  per  swine)  at  no  cost  to  this  project. 


5.2  Results  and  Discussion 

In  vitro  experiments 

Over  the  pressure  range  of  0  -  50  mmHg  (sim  ulating  the  IFP  in  breast  turn  ors)  OptiCell 
measurements  with  Definity  showed  an  inverse  linear  relationship  between  the  change  in 
subharmonic  amplitude  and  hydros  tatic  pressure  (r  2  =  0.79-0.99,  p  <  0.01).  This  is 
consistent  with  previous  results  reported  by  our  group  [13]  and  these  efforts  represent  the 
partial  fulfillment  of  tasks  la,  Id  and  If  in  the  original  Statement  of  Work  (SOW  ).  An 
example  of  the  subhannonic  am  plitude  for  Sonazoid  at  0  mmHg  and  47  mmHg  c  an  be 
seen  in  Fig  ure  1  (trans  mitting  fireq  uency  7.5  MHz  and  acoustic  pres  sure  0.7  M  Pa). 
However,  the  decrease  in  subharmonic  amplitudes  recorded  for  pressure  changes  from  0 
to  50  mm  Hg  varied  from  13.7  to  22.9  dB  a  nd  16.4  to  21.8  dB  (depending  on  acoustic 
pressures)  for  6.7  and  10  MHz  transm  ission  frequencies,  respectively.  These  levels  of 
subharmonic  am  plitude  variation  were  m  arkedly  higher  than  our  previous  results 
obtained  over  a  much  larger  pressure  range  (up  to  200  mmHg)  [13,  18-19]. 

After  som  e  extensive  testing,  we  discovered  that  standing  waves  m  arkedly  influenced 
results  and  made  measurements  impossible  to  reproduce.  Hence,  our  use  of  the  Opt  iCell 
setup  had  to  be  abandoned  and  a  com  pletely  new  tank  (with  extra  acoustic  abso  rbers 
incorporated)  had  to  be  designed.  T  he  new  water  tank  is  c  urrently  being  built,  bu  t  this 
issue  has  delayed  the  project  further  by  an  other  approxim  ately  4  months.  Thus,  the 
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project  is  approxim  ately  9  to  12  months  be  hind  schedule  and  we,  therefore,  intend  to 
request  a  one  year  no  cost  extension. 


Figure  1.  Comparison  of  the  subharmonic  amplitude  (circle)  measured  with  Sonazoid  at 
a)  0  mmHg  and  b)  47  mmHg pressures. 


The  previously  developed  sim  ulation  m  odel  of  the  dynam  ics  of  an  encapsulated 
microbubble  contrast  agent  was  modified  to  include  nonlinear  extensions  of  the 
viscoelasticity  and  this  model  has  now  been  published  [20].  The  intent  was  to  better 

account  for  the  experim  entally  observed  changes  in  subharmonic  signal  am  plitudes  as  a 
function  of  hydrostatic  pressure  s.  The  m  odel  showed  that  the  determ  ining  parameter  of 
subharmonic  response  is  the  ratio  of  the  excitation  frequency  to  the  resonance  frequency. 
Changing  the  am  bient  pressure  changes  th  e  resonance  frequenc  y  and  thereby  the 
frequency  ratio.  For  different  acoustic  ex  citation  pressure  levels,  changing  ambient 
pressure  can  either  increase  or  decrease  the  subharm  onic  response  depending  on  this 
ratio.  For  som  e  range  of  param  eters,  the  variation  is  far  m  ore  complicated.  T  his 
behavior  is  clearly  at  odds  with  the  expe  rimental  observations  m  entioned  above.  These 
discrepancies  m  ay  be  due  to  encapsulation  hue  kling  or  strain  softening  in  describing 
nonlinear  oscillations,  specifically  the  subharmonic  response  [20].  This  investigation  has 
been  submitted  for  publication  [21].  This  constitutes  the  continuation  of  tasks  lb,  lc  and 
le. 

In  vivo  experiments 

Software  to  analyze  RF  data  from  the  Sonix  RP  scanner  has  been  im  proved  and  the  best 
method  to  extract  the  subhann  onic  signal  co  mponents  from  the  frequency  spectrum  has 
been  estab  lished.  This  software  is  appli  cable  to  RF  data  acquired  from  any  tisssues 
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(including  breast  turn  ors).  However,  due  to  the  ava  ilability  of  ca  nine  ca  rdiac  data 
(obtained  with  funds  from  other  sources,  but  processed  with  funds  form  this  grant),  our 
initial  work  focused  on  cardiac  SHAPE  pressure  estimates  [22].  Briefly,  the  unprocessed 
RF  data  for  each  accumulated  pulse  (Figure  2A  )  was  transformed  to  the  Fourier  do  main 
and  the  subharm  onic  signal  am  plitude  (at  ha  If  the  funda  mental  frequency  -  Figure  2B) 
was  extracted  as  the  average  signal  in  a  40  %  bandwidth  around  the  subharm  onic 

frequency  (i.e.,  here  1.25  MHz).  The  extracted  subhannonic  signals  from  all  pulses  were 
processed  u  sing  a  m  oving  averag  e  filter  to  e  liminate  nois  e  spikes.  The  range  of  the 
subharmonic  signal  (i.e.,  the  difference  be  tween  maximum  and  m  inimum  subharmonic 
amplitude)  was  com  pared  from  each  pulse  cont  our  (after  elim  mating  excessively  n  oisy 
pulses)  for  each  incident  acous  tic  pres  sure.  Since  the  pressu  re  track  ing  using 
microbubbles  is  an  incident  acoustic  pre  ssure  dependent  phenom  enon,  the  incident 
acoustic  pressure  with  m  aximum  s  table  subh  armonic  rang  e  was  then  selected  for  LV 
pressure  tracking  as  shown  in  Figure  2C.  Results  from  four  canines  indicate  th  at  an 
overall  resolution  on  the  order  of  0. 19  to  5.48  mmHg  can  be  obtained  for  diastolic  left 
ventricular  pressures  if  the  aortic  pulse  pr  essure  values  are  known;  otherwise  resolution 
on  the  order  of  0.64  to  8.98  mmHg  can  be  obtained  relative  to  an  average  calibration 
standard  [  23].  Moreover,  a  pate  nt  applic  ation  has  been  subm  itted  based  on  this 
development  effort.  Thus,  considerable  tim  e  has  been  spent  on  optim  izing  this  software 
and  this  effort  represents  the  conclusion  of  task  2b. 

A  unique  opportunity  to  test  the  Stryker  pressure  m  onitoring  system  and  provide  proof- 
of-concept  of  the  use  of  SHAPE  for  estim  ating  IFP  was  pursued  in  th  e  Sinclair  swine 
melanoma  model.  Measurements  were  only  obtained  in  three  swine,  because  of  technical 
difficulties  (and  only  at  10  MHz  in  one  of  those  animals,  due  to  time  constraints).  At  6.7 
MHz  neithe  r  of  the  2  turn  ors  studied  showed  a  statistica  lly  sign  ilicant  relation  ship 
between  pressure  and  subharm  onic  signals  (p  =  0.2).  Most  likely  this  was  due  to  the 
LI 4-5  probe  having  too  high  a  bandwidth  an  d  the  L9-4  linea  r  array  was,  theref  ore, 
selected  for  the  rat  xenograft  studies  of  task  3.  Conversely,  at  10  MHz  one  tumor  showed 
a  subharmonic  decrease  of  1 1.3  dB  (r  2  =  0.90,  p  <  0.01)  and  9.7  dB  (r  2  =  0.82,  p  <  0.01) 
over  15  mmHg  for  acoustic  output  powers  of  -4  and  -8  dB  respectively  (Figure  3).  The 
second  animal  showed  a  decrease  of  7.3  dB  (r2  =  0.89,  p  <0  .01)  and  13.2  dB  (r2  =  0.98,  p 
<  0.01)  over  53  mmHg  of  variatio  n  (at -4  and -8  dB,  respect  ively).  This  difference 
between  tissue  and  tumor  IFP  is  rather  large  and  may  be  caused  by  the  sensitivity  of  the 
Stryker  system  to  the  angle  between  the  needle  and  the  tissu  e  being  studied  (this  issue  is 
currently  b  eing  inves  tigated  furth  er).  The  last  m  elanoma  s  wine  had  a  pressure 
differential  of  33  mmHg  between  the  turn  or  and  normal  tissue  with  a  corresponding  6.7 
dB  decrease  in  subharmonic  signal  am  plitude  at  a  -8  dB  acoustic  power  (r  =  0.92,  p  <0 
.04).  In  conclusion,  in  vivo  proof  of  concept  for  SHAPE  as  a  noninvasive  monitor  of  IFP 
has  been  provided;  albeit  still  based  on  a  ve  ry  small  sample  size  and  with  relatively  large 
subject-to-subject  variation. 
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Figure  2:  Steps  involved  in  extracting  and  processing  the  subharmonic  signal  for 
SHAPE.  A  typical  signal  obtained  with  pulsed  Doppler  (A)  and  the  frequency  domain 
representation  (B)  of  the  signal  in  (A).  The  processed  subharmonic  signal  from  all  the 
pulses  (solid  line)  and  the  pressure  catheter  data  (dotted  line)  are  shown  in  (C).  Note  the 
inverse  relationship  between  the  subharmonic  signal  and  the  pressure  obtained  via  the 
pressure  catheter  (in  agreement  with  our  previously  published  in  vitro  results). 


6  KEY  RESEARCH  ACCOMPLISHMENTS 

•  Hydrostatic  pressure  is  inversely  relate  d  to  the  c  hange  in  subharmonic  amplitude 
of  Definity  microbubbles  in  vitro  (r2  =  0.79-0.99,  p  <  0.01). 

•  SHAPE  experiments  were  conducted  in  vitro  in  an  OptiCell  setup,  but  standing 
waves  m  ade  results  difficult  to  reprodu  ce  and  a  new,  im  proved  setup  has  been 
designed. 
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Figure  3:  In  vivo  change  in  subharmonics  as  a  function  of  IFP.  While  the  slopes  are 
almost  identical  in  2  swine,  there  are  clearly  large  animal  to  animal  variations. 


•  A  com  puter  m  odel  to  sim  ulate  the  beha  vior  of  m  icrobubbles  as  a  function  of 
pressure  has  been  further  developed. 

•  Software  for  processing  of  in  vivo  SHAPE  data  has  been  optimized. 

•  In  vivo  proof  of  concept  for  SHAPE  as  a  noninvasive  m  onitor  of  IFP  has  been 
provided. 

•  The  L9-4  linear  array  has  been  selected  for  the  rat  xenograft  studies. 


7  REPORTABLE  OUTCOMES 

Publications 

V.  G.  Halldorsdottir,  L.  M.  Leodore,  B.  Cavanaugh,  F.  Forsberg.  Initial  in  vitro  study  of 
US  pressure  measurements  for  monitoring  neoadjuvant  chemotherapy  of  breast  cancer. 
Prog.  RSNA,  pp.  338,  2009. 

F.  Forsberg,  J.  Dave,  V.  G.  Halldorsdottir,  M.  McDonald,  J.  B.  Liu  C.  Leung,  K.  Dickie. 
Noninvasive  cardiac  subhannonic  pressure  estimation  in  vivo.  J  Ultrasound  Med,  vol  29, 
no  3,  pp.  S13,  2010. 
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F.  Forsberg.  The  holy  grail:  ultrasound  microbubbles  for  therapy.  J  Ultrasound  Med,  vol 
29,  no  3,  pp.  S17,  2010. 

F.  Forsberg.  Subharmonic  imaging  and  pressure  estimation.  Med  Phvs,  vol  37,  no  6,  pp 
3397-3398,2010. 

S.  Paul,  K.  Sarkar,  F.Forsberg.  Strain-softening  elasticity  model  of  the  encapsulation  of 
an  ultrasound  contrast  microbubble.  Am  Phys  Soc,  Div  Fluid Dynam, 
http://meetings.aps.Org/link/BAPS.2010.DFD.LR.l,  2010. 

A.  Katiyar,  K.  Sarkar,  F.Forsberg.  Subharmonic  response  from  ultrasound  contrast 
microbubbles  for  noninvasive  blood  pressure  estimation.  Am  Phys  Soc,  Div  Fluid 
Dynam,  http://meetings.aps.org/link/BAPS.20 1 0.DFD.HL.  1 ,  2010. 

J.  Dave,  V.  G.  Halldorsdottir,  J.  Eisenbrey,  M.  McDonald,  J.  B.  Liu,  J.  Raichlen,  C. 
Leung,  K.  Dickie,  F.  Forsberg.  Subharmonic  signals  for  noninvasive  cardiac  pressure 
estimation  -  initial  in  vivo  experience.  Accepted  for  publication  in  JUM,  2011 

*V.  G.  Halldorsdottir,  J.  K.  Dave,  J.  Eisenbrey,  P.  Machado,  J.  B.  Liu,  D.  A.  Merton,  B. 
C.  Cavanaugh,  F.  Forsberg.  Subharmonic  aided  pressure  estimation  for  monitoring 
interstitial  fluid  pressure  in  tumors:  in  vitro  and  in  vivo  proof  of  concept.  Accepted  for 
publication  in  JUM,  2011 

S.  Paul,  A.  Katiyar,  K.  Sarkar,  D.  Chatterjee,  W.  T.  Shi,  F.  Forsberg.  Material 
characterization  of  the  encapsulation  of  an  ultrasound  contrast  microbubble  and  its 
subharmonic  response:  strain-softening  interfacial  elasticity  model.  J  Acoust  Soc  Am. , 
vol.  127,  no.  6,  pp.  3846-3857,  2010. 

A.  Katiyar,  K.  Sarkar,  F.  Forsberg.  Modeling  subharmonic  response  from  contrast 
microbubbles  as  a  function  of  ambient  static  pressure.  Submitted  to  J  Acoust  Soc  Am, 
2010. 

V.  G.  Halldorsdottir,  J.  K.  Dave,  L.  M.  Leodore,  J.  R.  Eisenbrey,  S.  Park,  A.  L.  Hall,  Kai 
Thomenius,  F.  Forsberg.  Subharmonic  contrast  microbubble  signals  for  noninvasive 
pressure  estimation  under  static  and  dynamic  flow  conditions.  Submitted  to  Ultrasonic 
Imaging,  2010. 

J.  K.  Dave,  V.  G.  Halldorsdottir,  J.  R.  Eisenbrey,  J.  S.  Raichlen,  J.  B.  Liu,  M.  E. 
McDonald,  K.  Dickie,  C.  Leung,  F.  Forsberg.  Noninvasive  pressure  estimation  of  the  left 
ventricular  diastolic  pressures  using  subhannonic  emissions  from  microbubbles  -  an  in 
vivo  pilot  study.  Submitted  to  Circulation,  2010. 


Presentations 

October  8,  2009  QED  Program,  University  City  Science  Center,  Philadelphia, 

PA,  USA. 
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•  Non-invasive  pressure  estimation  using  ultrasound. 

October  22-23,  2009  24th  Annual  Advances  in  Contrast  Ultrasound  &  ICUS  Bubble 

Course  2009,  Chicago,  IL,  USA. 

•  In  vivo  subhannonic  pressure  estimation. 

The  95th  Scientific  Assembly  and  Annual  Meeting  of  the 
Radiological  Society  of  North  America,  Chicago,  IL,  USA. 

•  Initial  in  vitro  study  of  US  pressure  measurements  for 
monitoring  neoadjuvant  chemotherapy  of  breast  cancer. 

March  25  -  28,  2010  The  55th  Annual  Convention  of  the  American  Institute  of 

Ultrasound  in  Medicine,  San  Diego,  CA,  USA. 

•  Noninvasive  cardiac  subharmonic  pressure  estimation  in 

vivo. 

•  The  holy  grail:  ultrasound  microbubbles  for  therapy 

May  7,  2010  Spring  Symposium,  the  Delaware  Valley  Chapter  of  the 

American  Association  of  Physicists  in  Medicine,  Philadelphia, 
PA. 

•  Recent  advances  in  ultrasound  for  image  guided  therapy. 

July  18  -  22,  2010  52nd  Annual  Meeting  of  the  American  Association  of 

Physicists  in  Medicine,  Philadelphia,  PA,  USA. 

•  Subharmonic  imaging  and  pressure  estimation. 

September  30  -  25th  Annual  Advances  in  Contrast  Ultrasound  &  ICUS  Bubble 

October  1,  2010  Course  2010,  Chicago,  IL,  USA. 

•  In  vivo  cardiac  subharmonic  pressure  estimation. 

November  21-23,  2010  63rd  Annual  Meeting  of  the  American  Physical  Society, 

Division  of  Fluid  Dynamics  Long  Beach,  CA,  USA. 

•  Strain-softening  elasticity  model  of  the  encapsulation  of  an 

ultrasound  contrast  microbubble. 

•  Subharmonic  response  from  ultrasound  contrast 

microbubbles  for  noninvasive  blood  pressure  estimation. 

V.  G.  Halldorsdottir  was  selected  as  a  finalist  for  the  AIUM  2011  Young  Investigator 
Award  based  on  her  work  in  the  abstract  labled  *  above.  The  competition  will  be  judged  at 
the  upcoming  Annual  Conference  of  the  AIUM  (in  April  2011  in  New  York  City,  NY). 

8  CONCLUSIONS 

Definity  showed  an  inverse  linear  rela  tionship  between  the  change  in  subharmonic 
amplitude  and  hydrostatic  pressure  (r  2  =  0.79  -  0.99,  p  <  0.01)  over  the  pressure  range 
associated  with  breas  t  tumors  (0  -  50  mmHg)  when  m  easured  in  vitro  in  the  OptiCell 


November  29  - 

December  4,  2009 
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setup.  However,  we  discovered  that  standi  ng  waves  m  arkedly  influenced  results  and 
made  measurements  impossible  to  reproduce.  Hence,  our  use  of  the  OptiCell  setup  had 
to  be  abandoned  and  a  com  pletely  new  tank  (w  ith  extra  acoustic  abso  rbers  incorporated) 
have  been  designed  (currently  under  construction). 

Our  attem  pts  to  design  a  realistic  sim  ulation  model  acco  unting  for  the  experim  ental 
results  have  been  m  ixed  and  further  work  is  ongoing.  Software  for  analyzing  RF  data 
from  the  Sonix  RP  scanner  to  produce  SHAPE  pressure  estimates  has  been  successfully 
optimized  and  an  initial  publication  submitted  [23]. 

Finally,  in  vivo  proof  of  concept  for  SHAPE  as  a  noninvasive  monitor  of  IFP  (r2  >  0.81,  p 
<  0.01)  has  been  provided  in  a  swine  melanoma  model  (at  no  additional  expense  to  this 
grant).  This  work  was  selected  for  th  e  final  of  the  AIUM  2011  Young  Investigator 
Award  [24], 

In  summary,  task  1  has  been  partially  com  pleted  while  task  2  is  ongoing.  However,  due 
to  the  delay  caused  by  the  in  vitro  experiments  the  project  is  approxim  ately  9  to  12 
months  behind  schedule  and  we,  therefore,  intend  to  request  a  one  year  no  cost  extension. 
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Appendix  I 


The  Statement  of  Work  from  the  original  proposal: 

Objective  1 

Task  1:  Computer  modeling  and  in  vitro  experiments  (months  1-6) 

a.  Construct  an  in  vitro  experimental  pulse-echo  system  for  investigating  the  effect  of 
hydrostatic  pressure  variations  on  c  ontrast  m  icrobubbles  and  m  easuring  the 
resulting  changes  in  backscattering  (Month  1). 

b.  Design  and  m  odify  num  erical  co  des  for  a  theore  tical  mode  1  describing  the 
dynamics  of  contrast  microbubbles  under  different  pressure  conditions  (Months  1  - 
3). 

c.  Calculate  the  behavior  of  individual  contrast  microbubble  and  the  collective 
behavior  of  contrast  microbubble  populations  (Months  3  -  6). 

d.  Measure  changes  in  backscattered  fundamental,  second  and  subharmonic  signals 
for  different  contrast  agents  as  a  function  of  pressure  (Months  2  -  6). 

e.  Predict  optimal  contrast  agents  for  SHAPE  according  to  the  numerical  simulations 
(Month  6). 

f.  Select  optimal  contrast  agent(s)  for  SHI  and  SHAPE.  The  s  election  will  mainly  be 
based  on  experimental  measurements  (Month  6). 

Objectives  2-3 

Task  2:  Design  and  implementation  of  SHAPE  on  a  commercial  US  scanner  (months 
7-12) 

a.  Optimize  SHI  and  SHAPE,  based  on  in  vitro  measurements  and  simulations  using 
the  actual  parameters  of  the  designated  transducers  (Months  7-8). 

b.  Modify  a  state-of-the-art  US  imaging  system  (the  Sonix  RP)  to  incorporate  the  SHI 
contrast  imaging  modality  and  to  perfonn  SHAPE  (Months  8  -  10). 

c.  Evaluate  the  new  im  aging  modality  and  SHAPE  in  an  in  vitro  phantom  using  the 
modified  US  scanner  (Months  11  -  12). 

d.  Prepare  regulatory  review  and  obtain  approval  for  animal  studies  (Months  9  -  12). 

Objectives  3-5 

Task  3:  Animal  experiments,  data  collection  and  analysis  (months  13-  36) 

a.  Create  and  grow  breast  tumors  by  implanting  one  of  three  human  breast  cancer  cell 
lines  (SKBR3,  BT474  or  MCF-7)  into  the  mammary  fat  pad  of  athymic,  nude  rats 
(Months  13  -  34). 

b.  Calibrate  in  vivo  SHAPE  results  based  on  IFP  measurements  obtained  with  the 
intra-compartmental  pressure  monitor  in  21  nude  rats.  Three  groups  (one  per  cell 
line)  of  7  rats  with  breast  tumors  implanted  will  be  studied  (months  14  -  16). 
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c.  Produce  and  evaluate  the  ability  of  SHI  to  depict  normal  vascularity  as  well  as 
breast  tumor  angiogenesis  in  human  xenografts  implanted  in  nude  rats  compared  to 
CD31  stained  specimens  (Months  17  -  34). 

d.  Validate  the  clinical  potential  of  SHAPE  as  a  therapy  monitoring  tool  by  studying 
180  human  xenograft  breast  tumors  in  nude  rats  (42  normal  rats  and  138  after 
administration  of  a  chemotherapy  agent  paclitaxel)  and  comparing  results  to  intra- 
compartmental  pressure  measurements  (months  17  -  34). 

e.  Perform  statistical  analyses  and  write  final  report  (months  34  -  36). 
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